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Abstract 


An  escalation  method  for  performing  partial  fraction  expansions  is 
presented  for  the  case  that  the  complete  list  of  zeros  of  the  denominator 
of  the  proper  rational  function  is  known.  Expressions  for  the  number  of 
divisions  and  multiplications  required  are  developed.  The  new  method  re¬ 
quires  fewer  such  arithmetic  operations  than  does  the  method  of  Henrici. 

A  numerical  example  is  provided. 
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Introduction.  At  times  It  is  desired  to  express  a  rational  function  in  terms 
of  partial  fractions.  After  completely  factoring  the  denominator  polynomial 
into  linear  factors  it  is  conceptually  easy  to  perform  the  expansion,  but  more 
efficient  methods  of  carrying  out  the  calculation  are  always  welcome.  Particu¬ 
lar  attention  is  directed  to  the  situation  where  the  linear  factors  of  the  de¬ 
nominator  occur  repeated.  The  partial  fraction  expansion  technique  that  is 

presented  here  is  efficient  relative  to  the  number  of  arithmetic 
operations  required.  It  compares  favorably  in  this  regard  with  the  method  of 
Henricl  [2]. 

The  Escalation  Process.  Consider  the  proper  rational  function 


- u 

Q(x)(x-£a)  <x-£fe) 

where  £  ¥  £.,  A  and  B  are  integers  greater  than  zero,  and  P(x)  and  Q(x )  are 

3  D 

polynomials  for  which  neither  £a  nor  £^  are  zeros.  In  terms  of  partial  fractions 

♦  (X)  »  Y  — ^ — l  +  I  +  X(x>  (2) 

i-1  (x-cr  i-1  (x-c.) 

a  b 

Here  and  in  other  parts  of  this  paper  a  summation  is  taken  to  vanish  if  its 
upper  limit  is  less  than  its  lower  limit.  A  related  function  $(x)  has  the 
definition  and  partial  fraction  expansion  given  by 


a  i  /  \  A  C  ,  B 

♦  <*>  -  -  l  -Si-i  +  l  - 

'  V  i-1  (x-£  ) 1  i-1  (x-£t> 


In  going  from  $(*)  to  $(x)  the  power  of  (x-£a>  in  the  denominator  was  Increased 
by  unity.  For  this  reason  those  A  partial  fraction  coefficients  represented 
by  £  .  are  called  the  native  coefficients  of  (3).  All  other  partial  fraction 


1 


coefficients  of  (3)  (the  8^  and  those  contained  In  x(x))  are  called  alien  co¬ 
efficients.  What  follows  will  explore  and  exploit  the  conjecture  that  the  partial 
fraction  coefficients  (both  native  and  alien)  of  (3)  may  be  computed  from  their 
counterparts  in  (2). 

First  consider  the  equation 

(x-5.)A"1*(x)  -  (x-C  )A|(x)  (4) 

a  a 

and  the  (A-2)  equations  obtained  through  repeated  differentiation  with  respect  to 
x  (if  A  ■  1,  cr  2,  no  differentiation  is  indicated).  When  x  is  set  equal  to  £ 

a 

one  finds  that 


ai 


a,  i-1 


i  -  2,  3,  *•*,  A 


(5) 


No  information  is  gained  concerning  C 


al 


Next  consider 


(x-£b)B<Kx)  ■  (x-Cb)B(x-Ca)$(x)  (6) 

and  the  (B-l)  equations  obtained  through  successive  differentiation.  Upon 
setting  x-£b  one  finds  that 
ft  _  CbB 


C. — £ 
b  a 


and 


A 

»  C  -C 
£  bi  b.i+1 

T»i 

b  a 


B-l,  B-2, 


,2,1 


(7) 


(8) 


A  summary  is  now  presented  which  makes  use  of  signal  flow  graphs.  The 
second  row  of  coefficients  in  Figure  1  are  the  native  coefficients  of  (3)  while 
the  first  row  are  the  corresponding  coefficients  of  (2).  In  Figure  2  the  an¬ 
alogous  signal  flow  graph  for  the  computation  of  a  set  of  alien  coefficients  is 
presented.  The  procedure  must  be  repeated  for  each  different  set  of  alien  co¬ 
efficients.  For  brevity  take 
1 


f  - 


V«a 


(9) 


2 
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Application  of  Escalation.  Consider  the  proper  rational  function  given  by 


F(s)  = 


N(s) 

D(s) 


(10) 


where 


n-1 

N(s)  =  l  b  s1 

J  * 


(ID 


i=0 


and 


**  • 
D(s)  =  l  as1 
i=0 


(12) 


In  (12)  a  ^ 0  ,  but  no  such  analogous  restriction  is  implied  in  (11),  except 
n 

that  not  all  the  b.  are  zero.  The  distinct  zeros  of  D(s)  are  o  ,  o  ,  •••,  o 
i  1  2  q 

which  occur  with  multiplicities  of  M^,  M2,  ••*,  M  respectively.  The  complete 
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Figure  1.  Formation  of  native  coefficients. 


Figure  2.  Formation  of  alien  coefficients. 


list  of  zeros  of  D(s)  is  s^,  82*  •••,  s^  where  it  is  convenient,  but  not 

necessary,  to  arrange  the  list  such  that  the  first  items  all  equal  a^,  the 

next  M  items  all  equal  o„,  and  so  forth  until  the  final  H  items  all  equal  o 

2  2  <1  q 

In  factored  form  D(s)  becomes 


“(•)  -  *„  jSi  (—j)  -  *„  jlx  J 


(13) 


The  numerator  may  be  written  as 

n-1  i 

N(8)  “  e0  +  Bi  jSl  is'8j)  tt*> 

where  the  0  coefficients  are  found  by  Horner's  scheme.  The  details  of  this  cal¬ 
culation  will  be  illustrated  later.  Upon  defining 

Gn(s)  "  anF<s)  (15) 


one  obtains 


Vs> 


so  +  X  Bi  jil 


n  (s-s.) 


This  may  be  generalized  to 

Bo +  6i  jii<s-8j) 


G  (s) 
r 


n  (s-s.) 
j=l  j 


(16) 


r-1,  2, 


(17) 


Algebraic  manipulation  then  reveals 

Gr-l(s)  Br-1 
G  (s)  -  7-—-  .  +  -7--  1  , 

r  (s-s,)  (s— s  )  * 


r-1,  2, 


(18) 


provided  that  one  defines 


G0<»>  -  0  (19) 

All  of  the  Gr(s)  defined  by  (17)  are  proper  rational  functions  and  hence  have 
partial  fraction  expansions. 


The  immediate  aim  is  to  find  the  partial  fraction  expansion  of  Gr(s)  from 
the  corresponding  expansion  of  Gr_^(s).  If  this  can  be  accomplished,  it  can  be 
repeated  n  times  for  r=l,  2,  ...»  n,  ultimately  yielding  the  expansion  for  G^Cs) . 
Then  owing  to  (15)  the  expansion  of  F(s)  is  easily  found. 


The  process  of  going  from  the  partial  fraction  expansion  of  G^  ^(s)  to 
that  of  Gr(s)  may  be  broken  into  steps.  First  define 

A  0  (3) 

°r(s)--(fcrr  •  ,*1> 2-  — . »  <20) 

and  note  that  if  the  expansion  of  Gf  ^(s)  is  known  that,  with  the  exception  of 

A 

the  first  native  coefficient,  all  of  the  partial  fraction  coefficients  of  Gr(s) 
may  be  found  by  invoking  the  processes  shown  in  Figures  1  and  2.  Then,  ac- 


1 

cording  to  (18),  ~ is  added  to  the  foregoing  result.  This  addition  causes 

\S**8  / 


the  alteration  of  only  one  coefficient  in  Gr(s)  as  compared  to  G^/s).  That  altered 
coefficient  is  the  same  initial  native  coefficient  for  which  the  method  of  esca¬ 


lation  ^heds  no  light.  Thus  the  escalation  policy  as  given  in  Figures  1  and 


2  is  sufficient  to  transform  the  coefficients  of  Grl(s)  into  all  but  one  of  the 
coefficients  of  Gr(s).  The  remaining  initial  native  coefficient  of  G^(s)  may  be 
found  by  using  a  theorem  of  Hazony  [1]. 


Hazony's  theorem  states  that  the  sum  of  the  residues  of  any  rational  function 


whose  denominator  degree  exceeds  its  numerator  degree  by  two  or  more  is  zero.  It 


has  been  noted  that  G ^(s)  is  a  proper  rational  function.  This  means  that  the 

A 

denominator  degree  exceeds  the  numerator  degree  by  at  least  unity.  Since  Gr(s) 

A 

is  formed  according  to  (20)  it  is  clear  that  Gr(s)  fulfills  the  requirements  of 

A 

Hazony’s  theorem.  Thus  the  sum  of  the  residues  of  G^(s)  is  zero,  and  owing  to 
(18),  the  residue  sum  of  Gr(s)  is  Since  the  unknown  initial  native  coef¬ 

ficient  of  Gr(s)  is  a  residue,  its  value  may  be  found  by  noting  that  those 
partial  fraction  coefficients  of  Gr(s)  which  may  be  identified  as  residues 


Illustration  of  the  Method.  For  a  more  concrete  presentation,  attention  is 
directed  to  the  finding  of  the  partial  fraction  expansion  of 


F(s) 


l+2s+3s2+4s3+5s^+6s3 


24-104s+182s2-164s3+80s4-20s5+2s6 


(21) 


In  a  separate  (and  by  no  means  trivial)  calculation  it  may  be  found  that:  o^* 1 
M^*3;  02*2,  M2=2;  0^*3,  M^=l .  This  means  that  the  complete  list  of  the  zeros 
of  D(s)  is:  1,  1,  1,  2,  2,  3,  The  8  values  may  be  found  by  Horner's  scheme, 
which  involves  the  repeated  use  of  synthetic  division  using  the  complete  list 
of  zeros  of  D(s)  as  divisors.  The  initial  dividend  is  the  coefficients  of  N(s) 
and  the  remainders  are  the  8  coefficients.  The  process  in  abbreviated  form  is 
now  shown. 


ll 

[6  , 

5 

,  4 

,  3 

,  2 

ll 

[6  , 

11 

,  15 

,  18 

,  20] 

ll 

[6  , 

17 

,  32 

,  50] 

70 

ll 

(6  , 

23 

,  55] 

105 

ll 

[6  , 

35] 

125 

ll 

[6  ] 

47 

1] 

21 


It  may  be  verified  that 


N(s)  -  21  +  70(s-l)  +  105(s-l)  +  125(s-l) 

+47(s-l) 3(s-2)  +  6(s-l)3(s-2)2 
Next,  use  is  made  of  the  tableau  given  in  Figure  3. 


(22) 
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Figure  3.  Tableau  for  performing  escalation 


Each  row  of  the  main  body  of  the  tableau  contains  the  partial  fraction  co¬ 


efficients  of  G  (s)  defined  by 

q  Li 


Gj.  (s) 


I  l 


i=l  j=l  (s-o^' 


r-1,  2, 


n 


(23) 


where  the  L  are  the  appropriate  non-negative  integers  such  that  r  =  ][  L . , 

i=l  1 


where  none  of  the  exceeds  M^.  The  last  row  of  the  tableau  contains  the 


coefficients  for  G^(s) ,  which  are  easily  converted  into  the  coefficients  of 


F(s)  through  division  by  a 


n 


Those  regions  in  the  tableau  labled  "native  coefficients"  are  filled  in 
from  the  row  above  by  using  the  process  shown  in  Figure  1  and  by  demanding  that 
the  row  sum  of  those  coefficients  that  are  identified  as  residues  (they  bear 


check  marks)  sum  to  give  The  "alien  coefficients"  are  filled  in  from  the 


row  above  by  using  the  procedure  given  in  Figure  2.  The  calculation  progresses 
from  the  top  row  to  the  bottom  row. 

Omitting  the  details  of  the  calculation,  the  entries  in  the  tableau  are: 


21 

21 

70 

21 

70 

105 

70 

21 

105 

-196 

-91 

-21 

321 

125 

308 

112 

21 

-261 

321 

47 

1477 

8 

245 

4 

21 

2 

-60 

-321 

2005  , 

8  6 

i  of  the 

last  row 

coefficients  by  a  (=2) 
n 

gives  the  partial 

fraction  ex- 

of  F(s) . 

1477 

245 

21 

F(s)  = 

16 

"  8 

,  "  4 

(s-l) 

(s-l)2 

(s-l)3 

321 

2005 

+ 

_-30  + 

2  . 

16 

(24) 

(s-2)  + 

(s-2)‘ 

(8-3) 

aiiianftir  i 


Operations  count.  In  order  to  compare  the  efficiency  of  a  given  method  it 
is  useful  to  have  expressions  for  the  number  of  multiplications  and  divisions 
involved,  assuming  that  additions  and  subtractions  are  less  trouble¬ 
some.  Let  m  (<n)  be  the  degree  of  the  numerator.  In  order  to  find  the  g  co¬ 


efficients 


multiplications  are  required  when  Horner's  scheme  is  used. 


Referring  to  Figure  3  it  is  seen  that  the  native  coefficients  do  not  require  any 
multiplications  or  divisions.  This  assumes  that  the  multiplications  by  unity 
indicated  in  Figure  1  are  too  trivial  to  be  counted.  Each  alien  coefficient 
could  be  computed  with  only  one  division.  This  assumes  that  one  could  divide 

by  (o1-a  )  rather  than  multiplying  by  7 — - — r  as  indicated  in  Figure  2. 

1  n2  S  {<}l~a2) 

There  are  -j —  alien  coefficients  and  hence  that  many  divisions,  where 

S  -  M?  +  M2  +  * • •  +  M2  (25) 

1  Z  q 

Finally,  there  are  n  divisions  by  a^,  although  a  frequently  equals  unity.  The 
sum  of  all  the  above  operations  is 

n2  .  S  (m) (m+1) 

~2  +  n  ~  2  +  2 

The  last  term  may  range  from  zero  when  m  =  0,  up  to  when  m  =  n  -  1. 

Hence  the  operations  count  varies  from  a  low  value  of 
n2  S 

T  +  n-  2 

to  a  high  value  of 

2  n  S 

n  +  2  ~  2 

Henrici  [2]  gives  the  operations  count  for  his  method  as  being  less  than 
2 

2n  +  S. 

The  present  author  reanalyzed  the  method  and  under  the  assumption  that 


(n+1)  >  2M^,  i  ■  1,  2,  •••,  q,  found  that  the  low  value  (for  m=0)  is 
2n  +  n  -  2S 


9 


and  the  high  value  (for  m  *  n  -  1)  is 
3n2  +  -|n  -  |s 

Suppose  each  zero  were  to  occur  with  the  same  multiplicity  M,  then 
S  -  nM 

It  is  now  seen  that  when  n  is  large,  the  operations  count  for  escalation 
is  less  than  one  half  of  the  count  for  Henrici's  method. 
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